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In this letter we show that Quantum Vaccum Friction (QVF) should play an important role in 
neutron star evolution. Taking into account this effect we show that magnetars could be understood 
as a natural evolution of standard pulsars. For the Crab pulsar, of which the characteristic age is 
known, we present the first completely coherent time evolution for its period and braking index. 
For this pulsar we also give the predicted value of the current first derivative of the braking index, 
providing a very important test to confirm QVF. 
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One of the most fundamental Quantum ElectroDy- 
namics (QED) prediction concerns electric and magnetic 
properties of vacuum [l], Q. Following this idea, vac- 
uum can be regarded as a standard medium with its own 
electromagnetic properties. Nevertheless, effects related 
to these properties are extremely weak and can only be 
tested in the presence of a very high electromagnetic field. 
For that reason, pulsars are very appropriate systems to 
look for evidences for vacuum properties [3, 0] • Pulsars 
are fast-rotating neutrons stars, with a very high mag- 
netic dipole moment tilted with respect to their rota- 
tional axis. Typically, their mass is of the order of the 
solar mass, and their radius of the order of 10 km. But 
the most interesting feature of these stars is certainly 
their magnetic field which can exceed 10 8 T on their sur- 
face Pulsars are thus suitable laboratories to study 
magnetic properties of the surrounding vacuum. In par- 
ticular, a class of pulsars called magnetars [f| is charac- 
terized by a long spinning period and a magnetic field 
on the surface which is supposed to be even greater than 
the QED critical field 4 x 10 9 T @, Q. At such a field 
non linear properties of quantum vacuum cannot be ne- 
glected. In general, the surface magnetic field of pulsars 
is inferred by comparing the measured period derivative 
with the one obtained using the energy loss rate given by 
the classical formula for a radiating rotating magnetic 
dipole. We believe that inferring a field exceeding the 
QED critical one using a pure classical formula (which 
ne glec ts any QED effect) is questionable. Heyl et al. 
[jl. Il0| have shown that the energy loss rate of a radiat- 
ing rotating magnetic dipole, calculated by taking into 
account QED differs very little by the classical one. The 



use of a QED formula for the radiating dipole does not 
change in practise the inferred magnetic field on the sur- 
face of the pulsar. 

Very recently, we have shown that a new energy loss 
process for pulsars should be taken into account, namely 
Quantum Vacuum Friction (QVF) resulting from the 
interaction between the non-stationary magnetic dipole 
moment of the star and its induced quantum vacuum 
magnetic dipole moment (llj . In the case of magnetars, 
with very high magnetic fields, we have shown that the 
QVF is the dominating energy loss process and has im- 
portant consequences. In particular, the inferred value of 
the magnetic field no long exceeds the QED critical field 
for all known magnetars. 

In this letter we consider some astrophysical conse- 
quences of QVF. We show that QVF plays a crucial role 
in neutron star evolution. Taking into account this ef- 
fect, we show that magnetars can simply be understood 
as a natural evolution state of typical pulsars, similar to 
the Crab pulsar. For the Crab pulsar, of which the char- 
acteristic age is known, we present the first completely 
coherent time evolution for its period and braking index. 
Let us first summarize the main results concerning QVF 
in highly magnetized neutron stars. Considering a neu- 
tron star rotating in vacuum, we denote by R its radius 
and B the magnitude of the magnetic field at its surface. 
If P is the spinning period of the neutron star and 9 the 
inclination angle of its magnetic dipole moment with re- 
spect to the rotation axis of the star, one can show that 
for Bq < 10 10 T, the QVF energy loss rate E qv is given 
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where B c = m 2 e c 3 /eh ~ 4.4 10 9 T is the QED critical 
field, a the fine structure constant, /iq the vacuum per- 
meability and c the light velocity in vacuo. Total energy 
loss rate of the star E can thus be obtained by adding 
the classical dipolar radiative energy loss rate E r given 
by 
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It is important to stress that the spinning period depen- 
dence is not the same for the classical radiative process 
and for the QVF. The ratio of classical to QVF losses 
decreases like 1/P 2 for large P, so that QVF becomes 
more important for slowly rotating neutron stars. Let 
us now study how QVF drastically changes the evolution 
of a pulsar. The evolution of a pulsar in period P and 
slowdown P is usually described in a logarithmic P — P- 
diagram as shown in fig. [T] where we have plotted all of 
the currently known pulsars (~ 1700) for which P and 
P have been measured fl2j. On this diagram, we have 
also indicated lines of constant magnetic field in the case 
where only the classical radiative energy loss process is 
considered (dotted lines) together with the case for which 
QVF is taken into account (solid curves) . These lines are 
obtained by matching the loss in rotational energy E with 
E r + E qv . Since E is given by 
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where a neutron star moment of inertia of I — 10 45 g.cm 2 



is assumed, we get 
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where sin# is taken, as usual, equal to 1. The first term 
in eq. (j4|) corresponds to the classical contribution, which 
gives the asymptotic behaviour for small P. The second 
term is the correction due to QVF, which becomes im- 
portant for large P. In fig. [1] we clearly distinguish two 
different evolution areas : pulsars located in the upper 
right area of the figure (gray) have an energy loss pro- 
cess dominated by QVF, whereas pulsars in the lower left 
area have classical magnetic dipole emission as the main 
energy loss. Defining the borderline between both areas 
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Since P ^ 0, this equation implies that P = and thus, 
the borderline between both areas corresponds to pul- 
sars having a constant slowdown P. Assuming that the 
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FIG. 1: P — P-diagram for the currently known pulsars. Pul- 
sars located in the upper right area (gray) have an energy loss 
process dominated by QVF. The Crab pulsar is marked as a 
circle filled with a dot. 

pulsar's magnetic field to remain constant during its evo- 
lution, we can describe the evolution of the star follow- 
ing the line on which it is situated today. Let us first 
ignore QVF. Since the rotational period of a pulsar in- 
creases with time, the natural evolution of these stars 
corresponds to a fall towards the lower right area of the 
figure along a given straight dotted line, depending of 
the value of their magnetic field. In this case, pulsars 
located in the upper right area can exhibit surface mag- 
netic fields above the quantum critical field. Since these 
stars cannot be connected by a straight line of constant 
magnetic field to other pulsars with lower rotational pe- 
riod, we have to assume that these stars are very young, 
even if their rotational period is already high. For that 
reason, this class of pulsars, called magnetars, seems to 
be very different in nature from other standard pulsars. 
To understand the structure of these stars and the physi- 
cal process involved which should give rise to such a huge 
surface magnetic field remains an open question. If we 
take into account QVF, we see that all magnetars are 
now connected by a modified line of constant magnetic 
field to pulsars with lower rotational period. Remarkably 
enough, their surface magnetic field becomes lower than 
the critical quantum field. We can then conclude, as far 
as QVF exists, that magnetars could be understood as a 
natural evolution of standard pulsars. Such a very impor- 
tant astrophysical consequence of quantum vacuum effect 
shows that testing the validity of the QVF theory would 
be crucial in the future to understand pulsars physics. 
It is also important to stress that with magnetic fields 
below the critical quantum field, the current explanation 
for the emmission properties of Anomalous X-ray Pul- 
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sars (AXPs) and Soft Gamma-ray Repeaters (SGRs) can 
not be applied to our model [l3|, and so, this emmission 
remains an open question. Another interesting result of 
taking QVF into account in magnetized rotating neutron 
stars is a completely coherent time evolution of pulsars, 
as we will show in the following. Let us illustrate this 
point in detail with the Crab pulsar (the one marked as 
a circle filled with a dot in fig 111)- Situated in the center 
of the Crab Nebula which is the remnant of the Super 
Novae SN1054, the Crab pulsar is certainly one of the 



most well-known pulsars 14|. A very interesting feature 



of this pulsar is that historical records indicate its char- 
acteristic age to be ~ 955 yr. This indication, together 
with accurate measurements of its current spin period 
P, first and second time derivative of the spin period, re- 
spectively Pi and P2 , drastically constrains its rotational 
evolution dynamics. Rewriting eq. ((4]) as 
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where K r = E r x (P 4 /4tt 2 I) and K qv = E qv x (P 2 /Att 2 I), 
the age of a pulsar r is given by 
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where t denotes the current date and Pi the initial spin 
period at t = to- We obtain 
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If K qv <C K r eq. © becomes 
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From eq. (0 we have in this case K r ~ PP, so that 
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which is the expression usually employed to estimate the 
characteristic age of a pulsar. Inverting eq. © (resp. 
cq. (|lip) we obtain the spin period as a function of time 
when classical and QVF energy loss processes are taken 
into account (resp. when QVF is neglected). Both cases 
are presented in fig. [2] for the Crab pulsar. In this figure, 
t = corresponds to current time (i.e "today"). Since the 
age of this pulsar is known, we can estimate the initial 
spin period Pj. 

Using the classical expression we obtain Pj = 16ms, 
whereas we get Pj = 18 ms when QVF is taken into ac- 
count. Since the energy loss process of the Crab pul- 
sar is dominated by classical radiation (see fig. [IJ, we 
are not surprised that the global time evolution of this 
star does not drastically change when QVF effects are 
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FIG. 2: Spin period of the Crab pulsar as a function of time 
when QVF is taken into account (full line) and when QVF is 
neglected (dotted line), t = corresponds to current time (i.e 
" today" )- 

included. Nevertheless, the most important difference 
between both calculations is that the one taking into ac- 
count QVF is the only one which provides a spin period 
evolution completely coherent with current observations, 
in particular with the current value of the braking index, 
as we will discuss in the following. 
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FIG. 3: Braking index of the Crab pulsar as a function of time 
when QVF is taken into account (full line) and when QVF is 
neglected (dotted line), t = corresponds to current time. 

The braking index, n, is a fundamental dimensionless 
quantity given by |15| 
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where v is the spin frequency of the pulsar and v (resp. 
i>), denotes the first (resp. the second) derivative of v 
with respect to time. This parameter describes the rate 
at which a pulsar loses rotational energy. Introducing the 
spin period P = 1/V, we rewrite the braking index as 
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From this expression we see that the borderline between 
gray and white plotted in fig. [T] corresponds to pulsars 
having a braking index equal to 2 (P = 0). Pulsars 
located in the upper right area of the figure (gray) have 
a braking index lower than 2, whereas braking index of 
pulsars in the lower left area is higher. In the following 
we study the time evolution of n for t > 0. We denote 
by Pq the current spin period, and by P\ and P2, the 
current values of the first and the second derivative of P 
with respect to time. The current value of n is 
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At this stage it is important to stress that the constant 
value n = 3 is expected for pure classical radiation mod- 
els, whereas the current measured value is uq = 2.51 for 
the Crab pulsar. This discrepancy shows that pure classi- 
cal radiation models would not lead to a realistic picture 
for this pulsar's time evolution. Adding QVF effects, so- 
lution of eq. (O becomes, for t > 
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This expression allows us to calculate the first and the 
second derivative of P with respect to time, and to match 
P(0) with Pi and P(0) with P 2 , obtaining 



As far as we know, there are only two pulsars in the 
ATNF catalogue 12j for which the calculated value of 
n c (0) can be compared with our estimation given by 
eq. (|2"TT) . These stars are the Crab pulsar and the 
B1509-58 pulsar. For the Crab pulsar we get n(0) ~ 
-1.8421 lO" 1 ^- 1 and n c (0) 
for the B1509-58 pulsar we get h(0) 
a value in perfect agreement with the calculated one 
n c (0) ~ -1.2500 10~ n s _1 . These results show that new 
measurements of h(0) are needed to improve statistics 
and thus provide a very important test to confirm our 
QVF predictions. 

In this letter we have shown that QVF should play 
an important role in neutron star evolution. Taking into 
account this effect, we have shown that the surface mag- 
netic field of magnetars becomes lower than the critical 
quantum field, so that magnetars could be simply un- 
derstood as a natural evolution of standard pulsars. For 
pulsars for which the braking index is known, QVF pro- 
vides a completely coherent spindown history, as illus- 
trated with the Crab pulsar. For this pulsar we have also 
given the predicted value of the current first derivative 
of the braking index, suggesting a very important test to 
confirm QVF. 
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Finally, the time evolution of the spin period is given by 
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Using eq. (fl"8"|) . we can calculate the braking index 
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The braking index time evolution for the Crab pulsar is 
presented in fig. ©. When QVF is taken into account, 
the braking index is a decreasing function with respect to 
time which is consistent with a global rotation slowdown. 
From eq. (l20l) we can predict the current value of the first 
derivative of n given by 
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